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Abstract 

We employ the relativistic constituent quark model to give a unified description of 
the leptonic and semileptonic decays of pseudoscalar mesons (tt, K, D, D s , B, B s ). 
The calculated leptonic decay constants and form factors are found to be in good 
agreement with available experimental data and other approaches. We reproduce 
the results of spin- flavor symmetry in the heavy quark limit. 



1 Introduction 



Semileptonic decays of pseudoscalar mesons allow to evaluate the elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, which are fundamental 
parameters of the Standard Model. The decay K — > ireu provides the most 
accurate determination of V us , the semileptonic decays of D and B mesons, 
D — > K(K*)lv, B — > D(D*)lu and B — > n(p)lu, can be used to determine 
\V CS \, \Vcb\ and \V u b\, respectively. The effects of strong interactions in these 
processes can be expressed in terms of form factors, which depend on q 2 , the 
squared momentum transferred to the leptonic pair. Information on the form 
factors are obtained by measuring the distributions of q 2 and decay angles. 

The decays of heavy D and B mesons are of particular interest due to the spin- 
flavor symmetry observed for infinite quark masses [1]. This symmetry allows 
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to reduce the number of form factors and express them in terms of the uni- 
versal Isgur-Wise function [2]. Also the scaling laws derived for some physical 
observables can be, in principle, tested experimentally. Since the Isgur-Wise 
function cannot be calculated from first principles, many models and non- 
perturbative approaches, which exhibit the heavy quark symmetry, have been 
employed to describe relevant phenomena. However, it was found out, that the 
finite mass corrections are very important, especially, in the charm sector. It 
appears that in some sense a step back should be done from using the heavy 
quark symmetry as a guide under model building to the straightforward cal- 
culations with full quark propagators. Then one has to check the consistency 
of the results with the spin-flavor symmetry in the heavy quark limit. 

In this paper we employ the relativistic constituent quark model (RCQM) 
[3] for the simultaneous description of both light and heavy flavored me- 
son leptonic and semileptonic decays. This model is based on the effective 
Lagrangian describing the coupling of mesons with their quark constituents, 
and the compositeness condition. The physical processes are described by the 
one-loop quark diagrams with free constituent propagators and meson-quark 
vertices related to the Bethe-Salpeter amplitudes. The masses of lower-lying 
pseudoscalar (PS) mesons should be less than the sum of quark constituent 
masses to provide the absence of imaginary parts corresponding to quark pro- 
duction. The adjustable parameters, the widths of Bethe-Salpeter amplitudes 
in momentum space and constituent quark masses, are determined from the 
best fit of available experimental data and some lattice determinations. We 
found that our results are in good agreement with experimental data and other 
approaches. Also we reproduce the results of spin-flavor symmetry for leptonic 
decay constants and semileptonic form factors in the heavy quark limit. 

The shapes of vertex functions and quark propagators should be found from 
the Bethe-Salpeter and Dyson-Schwinger equations, respectively. This is pro- 
vided by the Dyson-Schwinger Equation (DSE) [4] studies. A DSE-approach 
has been employed to provide a unified and uniformly accurate description of 
light- and heavy-meson observables [5,6]. 

A similar approach, based on the effective heavy meson Lagrangian, has been 
described in [7] in terms of a model based on meson-quark interactions, where 
mesonic transition amplitudes are represented by diagrams with heavy mesons 
attached to quark loops. The free propagator has been used for light quarks. 
However, the quark propagator obtained in the heavy quark limit has been 
employed for heavy quarks. 
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2 The model 



We employ an approach [3] based on the effective interaction Lagrangian which 
describes the transition of hadron into quarks. For example, the transition of 
the meson H into its constituents qi and q 2 is given by the Lagrangian 

C- mt (x) = g H H(x) Jdx 1 Jdx 2 § H (x;x 1 ,x 2 )q(x 1 )T H \ H q(x 2 ) . (1) 



Here, Xh and are the Gell-Mann and Dirac matrices, respectively, which 
provide the flavor and spin numbers of the meson H. The function <3># is 
related to the scalar part of Bethe-Salpeter amplitude. For instance, the sep- 
arable form x±, x 2 ) = S(x — (x± + x 2 )/2)f((xi — X2) 2 ) has been used in 
[3] for pions. 

The coupling constants gu is given by the so called compositeness condition 
proposed in [8] and extensively used in [9]. That condition means that the 
renormalization constant of the meson field is equal to zero: 



Z H = 1 



4tt 2 



W H (m 2 H ) = 0, 



(2) 



where U' H is the derivative of the meson mass operator defined by 

•d 4 k 



IW) 



An 2 i 



T H S 2 (Jt)T H S 1 (f(+ j) 



(3) 



The invariant amplitudes describing the leptonic H(p) — > Iv and semileptonic 
H{p) — > H'{p')lu decays are written down 



A(H(p) -> ev) 



V2 



V qq ,(eO^)M»(p) 



G 



A(H(p) - H'(p')eu) = ^=V qq ,(eO„v)M£ H ,(p,p' 



(4) 
(5) 



where Gf is the Fermi weak-decay constant, V qq i is the appropriate element 
of the CKM matrix. The matrix elements of the hadronic currents are given 
by 



^S 2 (^S 1 ^+ p 1 ) 



f HP » (6) 
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M* H , (p, p') = -^9h9h> J ^<pH(-k 2 )<Pw(-k 2 ) 

l 5 S 3 (k') 1 5 S 2 (k'+ f)O^S 1 (^+ f) 



xtr 



= U(q 2 )(p + pr + f4q 2 )(p-pr 



(7) 



(8) 



where 4>h{— k 2 ) is related to the BS-amplitude in momentum space, and 

1 



Si(j6) = 



rrii— $ 



(9) 



is the propagator of the constituent quark with mass mj. As discussed before, 
to avoid the appearance of imaginary parts in Eqs. (6) and (7), we assume 
that rriH < m qi + m q2 which is a reliable approximation for the lower-lying 
mesons considered here. 

To evaluate the integral in Eq. (7) 

WP,P0 = /^^(-A; 2 )tr{ 7 5 53(^)7 5 5 2 (^+/)7^i(^+^)}, (10) 



where J-{— k 2 ) = (pH(—k 2 ) ■ (pH'(—k 2 ), we need to calculate the following inte- 
grals: 



j(0,fl,Hl/,IM/S) 



d 4 k {l,k»,k»k v ,Wk u k 5 )F{-k 2 ) 

7T 2 i [m 2 — (k + p) 2 ] \m\ — (k + p') 2 } [m 2 — k 2 ] 



Using the Cauchy representation for the function JF(— A; 2 ) and then the stan- 
dard techniques of the Feynman a— parametrization one finds {J-'{z) = dT{z) / dz) 



oo 2 3 

n 1> — 1 



;i2) 



oo „ 3 

n v i— 1 



(13) 



oo 2 3 



(14) 
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X 



H^TTi**'>-^(iTi) V <*'>} 



00 j. 2 3 

gTPi + ^ + 



(15) 



+ 



where q = p — p', P a = a±p + a 2 p', D 3 = aia 3 p 2 + a 2 a 3 p' 2 + aia 2 q 2 , and 
Finally, Eq. (10) becomes 

iUM) = ( P + p'Yi + (p 2 ,p' 2 ,q 2 ) + (p-p'Yi-(p 2 ,p' 2 ,q 2 ) 

with 



n 1— 1 



a, 



x 



4 - 3(«i + a 2 ) 



1 + t 



+ 



(mi + m 2 )m 3 

-(«i + a 2 )(mim3 + ^2^3 - mim 2 ) 



(16) 



+«i P 2 + a 2 p 2 ) - P a 2 ( 2 (2 - (ai + a 2 ) ^) ] } 



The normalization condition is written in the form 

3<& 



4tt 2 



/ + (p 2 ,p 2 ,0) = l 



(17) 



with mi = m 2 = m. 



The integrals corresponding to the matrix element of the leptonic decay H(p) — > 
Zz/ and radiative decay of neutral meson H{p) —> 7(51) + 7(52) are calculated 
following the same procedure. We have 
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y(p) 



An 2 i 



:0(-A; 2 )tr 7 5 S 2 (f6)^(I - 7 5 ) f) = P^{p 2 ) 



oo 1 
v y 



da 



m 2 + (m 1 - m 2 ) 



at 
1 + t 



<t>M (18) 



K^( qi ,q 2 ) = J — 0(-A; 2 )tr{ 7 5 W- ^7^)7^^+ 
= ie> xvap q«q 2 3 K{p 2 ) 



oo 2 1 1— ai 

K(p 2 )=mJ <*t(— ) / da i / da 2 (-0'(^) 





(19) 



where = t[a:m 2 + (l — a)??! 2 — a:p 2 + a: 2 p 2 t/(l + t)] and = £[m 2 — aia^p 2 ] + 
ai« 2 p 2 t/(l + 1). 

The physical observables are expressed in terms of the structural integrals 
written in Eqs. (16), (18) and (19): 



9Pyy 



9p 



2 v / 2tt ; 



7T 



-K(m 2 P ), T(P - 77) = -a 2 m^ 77 , 



(20) 



fp = — 2 g P Y{m 2 P ), T{P^lv) = \V q 



|2 P 



8tt 



-mpm l 



m. 



m, 



,(21) 



/+(? ) = 7T I+( m p, m p>,(l ), 



4tt 2 



r(P^P7z/) = |^ 



991 1927r 3 mi 



dt\f + {t)\' 



(f+-f)(t_-f) 



(22) 

-i 3/2 



with t± = (mp±mpr) 2 (the extra factor 1/2 appears for n° in the final state). 



2.1 Heavy quark limit 



The leptonic heavy decay constants and semileptonic heavy to heavy form 
factors acquire a simple form in the heavy quark limit, i. e. when m\ = M — > 
oo, m 2 = M' -> oo and p 2 = (M + £) 2 , p /2 = (M' + E) 2 with E being a 
constant value. From Eq. (16) by replacing the variables a\ — > a±/M and 
a 2 — > a 2 /M', one obtains 
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M + M' 

2MM' 



1 1 



- • [ dtf— 1 —) [ daa [ dr(-T'(z) 
' J \l + tj J J V w 



M + M' 1 )dr " 



2MM 





m + x /u 



m + 



at 

T+t 







m 2 + z 



(23) 



where z = u- 2E^/u/W, W = 1 + 2r(l - r)(w - 1) and w = (M 2 + M' 2 - 
2MM'g 2 )/(2MM'). 

The normalization condition can be obtained from Eq. (23) by putting w — 1 
and M' = M. We have 



3 9 H AO) _ 1 .(o) _ 1 , 



2M^' /7V 



y du^ 



(24) 



where Zq = u — 2E< s /u. Then the leptonic decay constant and semileptonic 
form factors are written as 



''^•l/S/^-^»<<?W (25) 

I N J W J Hy 1 m 2 + z K J 

v oo 

It is readily seen that we reproduce the scaling law for both leptonic decay 
constants and form factors, and obtain the explicit expression for the Isgur- 
Wise function [1,2]. 



3 Results and discussion 



The expressions obtained in the previous section for physical observables are 
valid for any vertex function <Ph{— k 2 ). Here, we choose a Gaussian form 
(f)(—k 2 ) = exp{k 2 / A H } in Minkowski space. The magnitude of A H charac- 
terizes the size of the BS-amplitude and is an adjustable parameter in our 
approach. Thus, we have six A-parameters plus the four quark masses, all of 
which are fixed via the least-squares fit to the observables measured experi- 
mentally or taken from lattice simulations (see, Table 1). 

The fit yields the values of A-parameters and the constituent quark masses 
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Table 1 

Calculated values of a range of observables (5^77 in GeV -1 , leptonic decay constants 
in GeV, form factors and ratios are dimensionless). The "Obs." are extracted from 
Refs. [10-16]. The quantities used in fitting our parameters are marked by "*". 



Obs. 



Calc. 



Obs. 



Calc. 



5 , 7T77 

u 

fK 



0.274 
0.131 
0.160 



0.242 
0.131 
0.160 



/no) 



0.98 

0.74 ± 0.03 



/ 



D 



0.191±^| 0.191 



is. 

Jd 
fB 

£b_ 

Sb 



1.08(8) 1.08 



0.2061^ 0.206 
0.172l^ 0.172 



' + 

/£*(0) 
Br (A" -»• irlv 



If? 1.14(8) 



1.14 
0.196 



0.27 ± 0.11 
(4.82 ±0.06) • 10~ 2 
Bv(D^Klu) (6.8 ±0.8) • 10- 2 
Bt(B->DIv) (2.00 ±0.25) • 10" 2 
Bt(B^ttIu) (1.8 ± 0.6) ■ 10- 4 



0.98 
0.74 
0.73 
0.51 

4.4 • 10~ 2 
8.1 • 10~ 2 
2.3 • 10~ 2 
2.1 • 10- 4 



which are listed in Eqs. (27) and (28). 



A. 



A 



K 



A 



D 



A 



D a 



A B A Ba (in GeV) 



1.16 1.82 1.87 1.95 2.16 2.27 



(27) 



TO,, 



TO., 



m c nib (in GeV) 



0.235 0.333 1.67 5.06 



(28) 



The values of A are such that A mi < A m . if m 8 < to,-. This corresponds to 
the ordering law for sizes of bound states. The values of Ad = 1.87 GeV and 
Ab = 2.16 GeV are larger than those obtained in [6]: A^ = 1.41 GeV and 
Ab = 1.65 GeV. The mass of u-quark and the parameter A n are almost fixed 
from the decays n — > \iv and 7r° — > 77 with an accuracy of a few percent. 
The obtained value of the u-quark mass m u = 0.235 GeV is less than the 
constituent-light-quark mass typically employed in quark models for baryon 
physics (m u > to^/3 = 0.313 GeV). For instance, the value of m u = 0.420 
GeV was extracted from fitting nucleon observables within our approach [3]. 
The different choice of constituent quark masses is a common feature of quark 
models with free propagators due to the lack of confinement. However, we 
consider here the low-lying mesons that allows us to fix the constituent quark 
masses in a self-consistent manner. As mentioned above, the meson masses 
must be less than the sum of masses of their constituents. This gives the 
restrictions on the choice of the meson binding energies: E K = m K — m s < m u , 
E D = m D — m c < m u and E B = m B — m h < m u , which means that the 
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binding energy cannot be relatively large as compared with those obtained in 
[6]: E D = 0.58 GeV and E B = 0.74 GeV. 

Let us now consider the g 2 -behaviour of the form factors. We use the three- 
parameter function for the four /+ form factors 

f HH'( 2\ = fjfi) ( 2 q\ 

I+ {q) 1 - W/m 2 H ) - &i(</7™ff) 2 1 ) 



here b , b\ and /(0) are parameters to be fitted. We collect the fitted values 
in the following Table and report the g 2 -dependence in Fig. 1. 





K -> 7T 


D -> K 


B -> D 


B -> 7T 




0.28 


0.64 


0.77 


0.52 


bi 


0.057 


0.20 


0.19 


0.38 



For comparison, we plot, together with our results, the predictions of a vector 
dominance monopole model: 

^ - rS^ (31) 



with m 2 / , being a mass of lower-lying go'- vector meson. We choose mp* = 2.11 

qq s 

GeV for c -> s, m B * = 5.325 GeV for 6 -> u, m B; w m Bc = 6.4 GeV [17] for 

6 — > c transitions. The values of /+ 9 (0) are taken from the Table 1. Also we 
calculate the branching ratios of semileptonic decays by using widely accepted 
values of the CKM matrix elements [10]. 

Our result for the slope of the K i3 form factor 

f K7r '(0) 

K = < Jk^ = °- 023 > (32) 



is in good ag reement with experiment: A+ xp = 0.0286 ± 0.0022 [10] and VDM 
prediction: A^ DM = m^/m 2 K * = 0.025. This value is also consistent with Refs. 
[18] 

One can see that the agreement with experimental data and lattice results is 
very good, with the exception of the value of /+ u (0) which is found to be larger 
than the monopole extrapolation of a lattice simulation, QCD Sum Rules (cf. 
[19]) and some other quark models (see, for example, [20,21]). However, this 
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result is consistent with the value calculated from Refs. [6,22] and allows us 
to reproduce the experimental data for B — > nlu with quite good accuracy. 
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q 2 (GeV 2 ) q 2 (GeV 2 ) 

Fig. 1. The semileptonic K—>tt,D—>K,B—>D and B — > it form factors with, for 
comparison, a vector dominance, monopole model Eq. (31) and a lattice simulation 
[23]. Our results: continuous lines. Monopole: dotted lines. Lattice: data points. 
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